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Abstract 


An improved isospin dependent Boltzmann Langevin model, in which the inelastic channels and 


momentum dependent interactions are incorporated, is used to investigate the high-density behavior of nuclear 


symmetry energy. By taking several forms of nuclear symmetry energy, we calculate the time evolutions of neutron 


over proton ratio, z multiplicity and z-/z* ratio, and the kinetic energy and transverse momentum spectra of z-/z* 


ratio in the heavy ion collisions at 4004 MeV. It is found that the neutron over proton ratio and z-/z^ ratio are very 


sensitive to the nuclear symmetry energy, and the z is more sensitive to the nuclear symmetry energy than the z*. A 


supersoft symmetry energy results ina larger z/z* ratio. 
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1 Introduction 


Constraints on the nuclear symmetry energy Esym(p) 
have been an very interesting subject in the last few 
years. The Esym(p) which is an estimate of the energy 
cost to convert all protons to neutrons at the fixed 
density in the symmetric nuclear matter, is very 
important both in nuclei, such as the binding energy, 
and in neutron stars, such as the nature and stability of 
the phases within the star! ?l. The energy per nucleon 
in asymmetric nuclear matter can usually be expressed 
as 

E(p,6) = E(p, = 0) E, (90^ +06"), — (1) 


where ó-(p,-p,)/p and p,, p, and p are the 


neutron, proton and total nucleon densities, 
respectively. The first term in Eq.(1) is the symmetric 
contribution and the second term is the symmetry 
energy contribution. The higher-order terms O(ó*) in 


ó are commonly negligible". Generally, there are 


two forms of Ess(p) predicted by various theoretical 
models. The one is that the E 


sym 
the increasing nucleon density and the other is that 


(p) increases with 


Esym(p) increases up to the saturation density p, and 
then decreases with the increasing density. 

At supra-saturation densities ( p > p, ), the form 
of the Esm(p) is very controversial. Even the trend 
of the Esm(p) with 
constrained yet. A lot of probes have been proposed to 


increasing density is not 


constrain the Esym(p) at supra-saturation density, 
such as the neutron over proton ratio of squeezed-out 
nucleons"!, the neutron-proton differential flow", the 
z [n ratio”, the K°/K* ratio!) the x-/x* 


HI and so on. Among above probes, the 2 /z* 


ratio 
ratio as a promising one is motivated by the A(1232) 
resonance model'?!, in which a primordial relation 
between z /z' ratio and N/Z is predicted. That 
is 


m [z' -(5N? + NZ (SZ! * NZ) «(N/Z) (2) 
where the N and Z are the neutron and proton 
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numbers in the participant region of the reaction. The 
N/Z ls Esym (p) by the 
dynamical isospin fractionation?. Therefore, one can 


determined by the 


use the z /z* ratio to constrain the Exs (p). 

Using the z/z 
Esym(p), a larger uncertainty exists. By comparing 
with the FOPI data?! a very soft Esym(p) was 
suggested by the isospin-dependent Boltzmann 
Uehling Uhlenbeck model (IBUU04)" and improved 
isospin dependent Boltzmann Langevin model 
(ImIBL)P), and a hard Ess(p) was suggested by the 
improved isospin dependent quantum molecular 
dynamics model (ImIQMD)"*!. 

At sub-saturation densities ( p < p, ), the trend 


ratio as a probe of the 


of the Esym(p) with increasing density is almost 
constrained but the quantitative results are not 
obtained yet. There are also a lot of probes to constrain 
the Esym(p) in sub-saturation density region, such as 
diffusion ^ ^! the 


thickness. the single and double neutron over proton 
[14,17-20 [21] 


the isospin neutron skin 


ratios | the heavy fragment mass distribution 


and so on. The IBUU04 modell”! suggested an almost 


linear Esym(p) by analyzing the isospin diffusion data. 


Also the constrained molecular dynamics II model 
(CoMDII) suggested a linear one by calculating the 


[21] 


heavy fragment mass distribution" . However, the 


improved quantum molecular dynamics model 
(ImQMD)!! and the isospin dependent quantum 
molecular dynamics model (IQMD)""7! predicted a soft 
Esm(p) by analyzing the double neutron over proton 
ratio. Recently, statistical 


a quantum approach 


predicted the Esym(p) at low temperature and low 
density limit""!, 

As a kind of ensemble-averaged theory, the 
usual BUU model have been very successful in 
describing the various phenomena in heavy ion 
collisions. However, the fluctuations, which are very 
important in preequilibrium particle 
subthreshold 


emission, 


particle production and 


multifragmentation are not included in the BUU model. 


In the QMD-like model which is a kind of many-body 
theory, the reactions are simulated event by event. 
Though the fluctuations are realized in the QMD 
simulations, the correlations among the events are 
uncertain. Moreover, the QMD model fails to explain 


the subthreshold kaon production""l, 


In this paper, we will use the ImIBL model to 
investigate the Esm(p) using the neutron over proton 
ratio and z /z' ratio in the "Ca Ca collisions at 
4004 MeV. The paper is organized as follows. We 
introduce our model in Section 2. The calculated 
results are given in Section 3. Finally, we will give 
some conclusions. 


2 Theoretical model 


The ImIBL model used in this paper is an updated 
version of the Boltzmann-Langevin model(BL)"*l. 
According to the BL model, the fluctuating phase 
space single particle density fulfills the following 


equation!”"!: 


9 p / f (F, B,t)= 
(a Evry, rna 6) 


K(f 2k (. p.1). 


The left-hand side describes the Vlasov propagation 
determined by nuclear mean field U ( f ). K ( f ) is 
the collision term of the usual BUU form but is 
characterized by the fluctuating density f (F, p,t). 
OK (7, p,t) is the fluctuating collision term, which can 
be explained as a stochastic force acting on the 
7 (7,p,t), and is characterized by the following 
correlation function: 


(OK (F, Bist, OK (5. Bast, )) = 
C(B,. D;)ó(, —7,)6(t, =t), 


where the angle brackets denote a local average, 


(4) 


performed over fluctuating densities generated during 
a time interval ôt. The correlation function C(p,, p) 
can be reduced and expressed as in the weak-coupling 
limit and determined by the one-body properties of the 
locally averaged distribution as indicated in Ref.[25]. 
The numerical simulation method used in BL model is 


the projection method? ^?! 


which projects the 
fluctuations on a set of low order local multipole 
moments Qj, of order L with magnetic quantum 
numbers M of the momentum distribution. The 
fluctuations of these multipole moments can be 


characterized by a diffusion matrix and can be 
[25] 


> oy 


deduced from the correlation matrix C(p,p') as 
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with 
AQ = On, (À.)* Qi, (5;.)- Qu, (B;)- Quy (24). 


According to the Refs.[24-27], the stochastic 
component of the dynamics can be reduced to the 


following two equations: 


O,)(7,t+ 5t) =O, (7,t+6t 
l (6) 


+ JótC, (7,1)W,, 


) 
Oy (7,t+ ôt) = Q, (F,t+ 8t) + ôt, (7,0), 


where W, and W; are Gaussian random numbers. 
The Q) and Q, are the fluctuating quadrupole and 
octupole moments. The fluctuations are inserted into 
the momentum space through a scaling procedure as 
indicated in Ref.[25]. 


“cat“Ca 
E =100A MeV ^ ^ ^ 


10 20 30 40 50 60 70 80 90 100 
t/ fmc! 
Fig.1 (Color online) Time evolution of the ensemble-averaged 
quadrupole moment Q4 , Q,+0, and Q,*2o, of the 
momentum distribution in the ImIBL model (upper) and the 
associated variance o, (bottom) for central “°Ca+“8Ca 
collisions at 1004 MeV. 

In Fig.l, we show the time evolution of the 
ensemble-averaged total quadrupole moment Q20, 
Q,,+0, and Q,+20, of the 
distribution and the associated variance o,, from 
central “Cat+*8Ca collisions at 1004 MeV in the 
ImIBL and BUU models, where o» = (03) - (Qu) 
is the standard deviation function. For a Gaussian 


momentum 


distribution, Q,, o, and Q,,+2o0,, correspond to 


84.3 and 99.5 percent of the number of events 
respectively". From the upper panel of Fig.l, we can 
see that the Q,, value of the ImIBL model has a 
larger range of variation as compared with the BUU 
model. From the bottom panel of Fig.l, there is a 
bump in the early collision stage for BL model. But 
this phenomenon does not happen in the BUU 
simulations. This phenomenon suggests that there is a 
larger difference between the real quadrupole moment 
and averaged quadrupole moment, especially in the 
early collision stage. We think this phenomenon is 
important in the dynamical evolution process. 

100 


—— |mIBL HP 
---ImlQMD hard, ; 
e ImQMD £3 


PI Po 


Fig.2 (Color online) Symmetry energy as a function of the 
reduced density in this work (solid line), Ref.[8] (dashed line) 
and Ref.[14] (dotted line). 

In the Vlasov part of Eq.(3), we work with an 
isospin and momentum dependent single nucleon 
potential, which reads 


v (pp) en [ | e Eis (p)? 


0 Do 
OES ( ) 38? 
OR 90? +E — (7 
2 sa (P)P 23 (7) 
+U jinis 


the bulk parameters a, f, and y taken here are 
—390 MeV, 320 MeV and 1.14 for the soft EOS plus 
MDI as SM and -130 MeV, 59 MeV and 2.09 for the 
hard EOS plus MDI as HM"? The compressibilities 
K are 200 and 380 MeV for the SM and HM, 
respectively. The E» 


the Esm(p), taken here are the following two forms 


" 1 Ys 
Era (A) = 5 Com n (8) 


(p), which is the local part of 
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and 


2 5/3 
gt) +2) . (9) 
Po Po Po 


In Eq.(8), the y,=0.5, 1.0 and 2.0 correspond to the 
soft, linear and hard symmetry energy respectively. 
The parameters Cym, a, b and c are 33.4, 38.9, 
—18.9, and —3.8 MeV, respectively. Eq.(9) gives the 
supersoft symmetry energy. The curvature parameters 
Kym 799) (9 E, fap"). are 2752, 254, —62.9, 
—403.6 MeV for the hard, linear, soft, and supersoft 
symmetry energy, respectively. Shown in Fig.2 is the 


comparison of the Esm(p) for different cases from 
Eqs.(8) and (9), and also from the Refs.[8,14]. 
The momentum dependent potential can be 


expressed asl 


Uy, 715710 | 0.0005( 3, E A “fe. (10) 

0 
The inelastic channels and the parameterized 
cross sections of each channel are taken from Ref.[29]. 
We take the free space and in-medium cross sections 
for the elastic and inelastic channels, respectively. The 
in-medium effects are the medium correction of p 


30 
mass! 1, 


3 Results and discussion 


3.1 Neutron over proton ratio 


Figure 3 shows the time evolution of the free neutron 
over proton ratio and the corresponding central point 
density in central “*Ca+**Ca collisions at 4004 MeV. 
The calculations are performed by using the supersoft 
and hard symmetry energies, respectively. The solid 
and dashed lines denote the results in the angle 
interval 70° «0 «110? and in the rapidity interval 
—0.2« y? «0.2. The dotted and dotted-dashed lines 
show the results without any cuts. 

From Fig.3, we can see that the supersoft 
behavior of the symmetry energy leads to a low 
neutron emission in the high-density region. In 
low-density region, an inversion of this trend is 
Observed because of contributions to the nucleon 
emission in the expansion phase while fragments are 
forming. Our results are consistent with the one of 
Ref.[31]. 


—— supersoft '70°< 9«110 
- - -hard 70°< 9 «110 


—-—- Supersoft 


0 —_—- 
0 10 20 30 40 50 60 70 
t/ fmc? 


Fig.3 (Color online) (Upper) Time evolution of free neutron 
over proton ratio for central ^5Ca--"*Ca collisions at 4004 MeV 
in the central rapidity and in the emission angle interval as 
indicated. (Bottom) Time evolution of the corresponding central 
point density. 


3.2 Pion multiplicity 


In the left panel of Fig.4, we show the time evolution 
of z^ multiplicity in central "Ca-* Ca collisions at 
400A MeV. The calculations are performed by taking 
four different forms of the Esym(p). From the left 
panel of Fig.4 and the bottom panel of Fig.3, one can 
see that the z is mainly produced in the high density 
region. Therefore, the z production will take useful 
information of the high density phase space. Most 
interestingly, the produced z is more sensitive to 
the Esm(p) than the z'. In order to further shed 
light on this point, we show in the right panel of Fig.4 
the z^ multiplicity as a function of the curvature 
parameter K.m. From the right panel of Fig.5, we can 
see that the z yield strongly varies with different 
Esm(p) from soft to stiff case. This is because the 
z is mainly produced in neutron-neutron collisions. 
In the high density region a softer Esym(p) will result 
in a low neutron emission, and then there will be more 
neutron rich as compared with the case of a hard 
Esm(p). We can conclude that the z^ yield play an 
important role in investigating the Esym(p). 


3.3 m/m ratio 


To further understand the dependence of the Esym(p) 
on the z production, we show in Fig.5 the 
differential z /x* ratios as functions of the kinetic 
energy and transverse momentum in the ^ Cat" Ca 


collisions at 4004 MeV with the impact parameter 
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b —] fm. One can see that in the low energy and low 
transverse momentum regions the z /z^ ratio is very 
sensitive to the Ess(p) . Therefore, we do not need to 
measure the high energy and high transverse 
momentum pions to constrain the density dependence 
of the symmetry energy. This will be useful in 
designing the experimental facilities with an aim of 
investigating the density dependent symmetry energy. 
From Fig.5, it is found that a supersoft symmetry 
energy results in a larger a [n ratio, which is 
of Refs.[7,9] and 
inconsistent with the one of Refs.[8,10]. 


consistent with the results 


04 0.28. 


M(x) 
r 


Mr) 


-400 -200 0 200 400 


t| fmc? Kom | MeV 


Fig.4 (Color online) (Left) Time evolution of z^ 
multiplicity and (Right) the z^ multiplicity as a function of 
the K,, in the “*Ca+“*Ca collisions at 4004 MeV with the 


sym 


impact parameter b =1 fm. 
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Fig.5 (Color online) (Left) Kinetic energy spectra and (Right) 
transverse momentum spectra of the zr /z' ratio for the 
^*Ca-^*Ca collisions at 4004 MeV with the impact parameter 
b =1 fm. 


4 Conclusion 


In summary, the Esm(p) has been investigated by 
using the free neutron over proton ratio and the z 
production in the ^5Ca--5Ca collisions at 4004 MeV in 
the framework of ImIBL model. It is found that both 
the free neutron over proton ratio and the x /z^ ratio 
are sensitive to the Esym(p). A soft Esym(p) results in 


a low emission of neutron component in the high 
density region. The z yield is strongly dependent on 
the Esm(p) and the z^ yield is almost independent 
on the Exs(p). The m` / z' ratio in the low kinetic 
energy and low transverse momentum region is 
sufficient to constrain the Esym( p). The z | z* ratio 
becomes lager when the Esym(p) varies from the stiff 
to the soft. 

Meson production as a probe of the Esym(p) is 
still an open question. It is known that the threshold 
effect is important in the K meson production. For 
z meson, the threshold effect, which is not concluded 
in the present calculations, may also play an important 
role in the whole evolution stage. We will incorporate 
the threshold effect into the ImIBL model in the future 
work. 
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